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Abstract.  "The  uniqueness  theorem  of  this  paper  answers  an 
open  question  for  a  system  of  differential  equations  arising  in; 
a  certain  n-body  problem  of  classical  electrodynamics.  The 
essence  of  the  result  can  be  illustrated  using  the  scalar  proto¬ 
type  equation  x'  ■  g-^x)  +  62^  +  x)  with  x(0)  *  0.  The 
solution  of  the  latter  will  be  unique  provided  g-^  and  g2  are 
continuous  positive  functions  of  bounded  variation .r 
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The  theorem  proved  in  this  paper  presents  a  criterion 
weaker  than  a  Lipschitz  condition  which  assures  uniqueness  of 
solutions  of  a  system  of  ordinary  differential  equations.  It 
was  designed  to  resolve  an  open  question  in  classical  electro¬ 
dynamics  described  at  the  end  of  the  paper. 

Before  stating  the  theorem  let  us  illustrate  it  with  two 
scalar  examples  typifying  the  problems  we  had  in  mind.  These 
examples  are  easily  treated  with  the  theorem  which  follows. 

We  are  unaware  of  any  previous  uniqueness  theorem  which  would 
handle  them  or  the  electrodynamics  problem  of  Example  3. 

Example  1.  If  g^  and  g2  are  continuous  positive 
functions  of  bounded  variation  on  an  open  interval  containing 
0,  then  the  equation 

x*  =  g^x)  +  g2(t  +  x)  with  x(0)  *  0 
has  a  unique  solution  on  some  open  interval  containing  0. 


Example  2.  The  equation 

x’  =  (t  +  x5/3)1/3  for  t  ^  0  with  x(0)  =  0 
has  a  unique  solution. 


The  theorem  itself 
differential  equations 
(1)  x’ 

with  initial  conditions 


(2) 
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Let  S  be  a  subset  (not  necessarily  open)  of  Rn  x, 
and  let  f:  S  -*•  Rn.  Then,  given  (tQ,  xQ)  £  S,  a  solution  of 
Eqs.  (1)  and  (2)  is  defined  as  any  differentiable  function 
x  on  an  interval  J  such  that  (t,  x(t))  €  S  and 
x'  =  f(t,  x(t))  for  t  £  J,  while  tQ  €.  J  and  x(tQ)  =  xQ. 
(If  J  contains  either  of  its  endpoints,  x'(t)  is  a  one-sided 
derivative  there . ) 

The  norm  used  in  this  paper  for  a  vector  5  €  Rn  is 

lull  =  Uil. 


Theorem.  Let  f:  S  R  be  continuous  and  satisfy  the 
following  condition.  Each  point  in  S  has  an  open  neighbor¬ 
hood  U,  a  constant  K  >  0,  an  integer  m  >.  0,  and  functions 
hj  and  gj  for  j  *  1,  ...,  m  such  that 

m 

(3)  | |f(t,c)  -  f(t,n) | |  <  K| \K  -  n| |  +  K  E  Jgj  (hj  (t,£) )  -  gj  (hj  (t,n) )  | 


on  HA  S,  where  h^ :  U  -*•  R  is  continuously  differentiable 
with 

3h,(t,C)  n  3h.(t,£) 

(4)  — 1 ^ -  +  Z  ^  —  f^t,?)  t  0  on  UA  s 


and  each  gj  :  R  R  is  continuous  and  is  of  bounded  variation 
on  bounded  subintervals.  Then  Eqs.  (1)  and  (2)  with  any  point 
(tQ,  Xq)  £  S  have  at  most  one  solution  on  any  interval  J. 


Remarks.  The  theorem  of  course  does  not  guarantee  the 
existence  of  a  solution  on  a  nontrivial  interval  J.  Existence 
would  follow,  for  example,  if  S  were  open. 

To  treat  Example  1,  define  h^t,  £)  *  S  and  h2(t,  5) 

■  t  +  5.  For  Example  2,  let  h(t,  5)  *  t  +  and  g(£)  ■  5^^. 
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Proof  of  the  Theorem.  Suppose  there  were  two  different 
solutions,  x  and  y,  on  some  interval  J  =  [t^,  b)  where 
b  >  tQ.  (The  case  J  =  (b,  tQ]  is  handled  similarly.)  Let 
a  =  inf  (t  €  (tQ,  b)  :  x(t)  #  y(t)K 
Then  x(a)  =  y(a) . 

For  the  poifrt  (a,  x(a))  £  S  let  U,  K,  m,  hj ,  and  gj  be 
as  described  in  the  hypotheses  of  the  theorem.  Without  loss 
of  generality,  assume  that  for  each  j  the  expression  in  (4) 
is  positive  at  (a,  x(a)).  Then,  reducing  U  if  necessary, 
the  continuity  of  the  derivatives  of  hj  assures  that  there 
exist  positive  constants  p  and  M  such  that  for  j  3  1,  ...»  m 
3h,(t,€)  n  3h,(t,C) 

( 5 )  — “•  —  +  Z  — —  f ^ ( t , £ )  >_  p  on  UAS 

and 

(6)  |hj(t,€)-hj(t,n)|<M|U-n||  on  U. 

Choose  a  bounded  interval  [<*j ,  0j]  which  contains  hj(U/lS), 
reducing  U  if  necessary.  Then  gj  is  the  difference  of  two 
continuous  non-decreasing  functions  on  [<jj ,  6 j ] ,  and  each  of 
the  latter  can  be  extended  to  a  continuous  non-decreasing 
function  on  R  by  defining  it  to  be  constant  on  (-»,  a  j  ]  and 
constant  on  [gj,  ») .  Without  loss  of  generality,  we  shall  assume 
that  each  gj  is  itself  non-decreasing  on  R  and  that 
(7  )  gj (hj (a,  x(a)))  *  0. 

Define 

z(t)  *  f  ||x*(s)  -  y’(s)j|  ds  for  a  <  t  <  b . 

J  a 

Then  z(a)  *  0,  z'(a)  *  0,  z  and  z'  are  continuous, 
z’(t)  >  0,  and  | | x ( t )  -  y ( fc )  |  |  <  z(t)  on  [a,  b). 
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Choose  c  £  (a,  b]  sufficiently  small  so  that  (s,  x(s)) 
and  (s,  y(s))  remain  in  U  for  a  £  s  <  c.  Then,  from  (6), 
hj(s,  x ( s ) )  -  Mz(s)  <  hj(s,  y(s))  <  hj(s,  x(s))  +  Mz(s) 

and ,  from  ( 5 ) , 

fg-  hj  (s,  x ( s ) )  >  p 

for  a  s  <  c  and  j  =  1,  m. 

Thus  for  a  ^  t  <  c,  using  (3)  and  the  mono tonicity  of  each  g ^ , 

z(t)  <  K  f  ||  | x(s )  -  y (s )  |  |  +  E^|gj(hj(s,x(s)))  -gj(hj(s,y(s)))||ds 

<  K(t  -  a)z(t)  +  £  E  [g,(h,(s,x(s))  +  Mz(s)) 

P  j  =i J  a  i  " 

-  gj  (hj  (s,x(s) )  -Mz(s))]  3s"  hj(s,x(s))  ds 

m  /-hj  (t,x(t)  )+Mz(t) 

=  K(t  -  a)z(t)  +  ^  E  gj(u)  du 

^  1  hj  (t,x(t)  )-Mz(t) 

*#•  in  / 1 

-  -  E  I  [gj(hj(s,x(s) )+Mz( s) )  +  gj (hj (s ,x(s) )-Mz(s ) ) ]Mz' (s )ds . 
j  =1  ^  3. 

Choose  6^  >  0  such  that  for  each  J 

|gj(u)|  <  when  |u  -  hj(a,x(a))|  <  6^ 

Then  choose  6  £  (0,1/6K)  such  that  a  +  $  <  c  and, for  each  j, 
[hj(t,x(t))  -  hj(a,x(a))|  +  Mz(t)  <  when  a  <  t  <  a  +  d. 

Now  for  a  <  t  <  a  +  5  one  finds  z(t)  <  5z(t)/6.  This 
contradiction  completes  the  proof. 

The  motivation  for  this  paper  was  the  following  problem 
from  classical  electrodynamics. 
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Exaraple  3.  Consider  n  electrically  charged  point 
particles  moving  along  the  x-axis  at  distinct  positions, 
x-^t),  x2(t),  xn  ( t ) .  Assume  that  the  motion  of  particle 

j  depends  only  on  the  electromagnetic  fields  produced  by 
the  other  n-1  particles,  with  these  fields  traveling  to 
particle  j  at  the  speed  of  light,  c. 

The  required  fields  are  calculated  in  terms  of  the 
trajectories  of  the  other  particles  from  the  retarded 
Lienard-Wiechert  potentials;  and  they  are  substituted  into 
the  Lorentz  force  law  for  particle  j.  Introducing  v1  =  x^/c 
for  the  velocity  of  particle  i  as  a  multiple  of  c,  one 
obtains  a  system  of  delay  differential  equations  with  state- 


dependent  delays: 


(1  -  v])3'2 


l  KU  +  vl(t  -  rlJ>  _ 

W  vli  -  vi(t  '  riJ>  ’ 


where  each  K^j  is  a  constant,  =  sgn  [Xj(0)  -  x^(0)]. 


and  where  r^  >  0  satisfies 


r.  .  'i  -  V*  -  ru> 

iJ  0ij  -  vi(t  -  ru> 


for  i  *  j. 


In  these  equations  v^  and  r^j  without  an  argument  stand 
for  v^(t)  and  r^Ct). 

2 

In  order  to  solve  the  system  of  n  equations  represented 
by  (8)  and  (9)  when  t  >  0,  one  should  know  not  only 


v j ( 0 )  and 


for  all  J  and  all  i  j*  j. 


but  also  the  values  of  v^t)  for  t  <  0,  i  ■  1,  ...,n. 
Now  consideration  of  the  problem  in  three-dimensional 


motion  has  led  to  the  conclusion  that  accelerations  should  not 
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be  assumed  continuous,  but  only  integrable  [2],  Thus  it  seems 
reasonable  even  in  the  case  of  one-dimensional  motion  to  assume 


that  the  given  past  history  of  v^,  say 

(11)  VjCt)  =  gt(t)  for  t  <  0  (i  =  1,  n) 

is  merely  absolutely  continuous — not,  in  general,  locally 

Lipschitzian. 

Substituting  (11)  into  the  right  hand  sides  of  Eqs .  (8) 
and  (9)  one  gets  a  system  of  ordinary  differential  equations 
which  satisfies  the  uniqueness  criterion  of  the  present  paper. 


Thus  a  unique  solution  exists  at  least  as  long  as  each 


t  -  r^j(t)  ^  0  and  each 
of  the  solution  would  use  a 


|vj(t)|  <  1.  (Further  extension 
"method-of-steps"  argument  which 


is  not  relevant  to  this  paper.) 


The  above  uniqueness  problem  was  solved  earlier  for  the 
case  of  two  particles  in  one-dimensional  motion  [1].  But  the 
method  used  did  not  seem  to  extend  to  the  n-body  problem. 
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